A0-A102  905 

UNCLASSIFIED 


AIR  FORCE  GEOPHYSICS  LAB  HANSCOM  AFB  MA  F/g  l^/i 

APPROXIMATE  ANALYTIC  SOLUTIONS  FOR  THE  PRIMARY  AURORAL  ELECTRON— ETC (U> 
MAR  01  d  R  JASPERSE»  D  J  STRICKLAND 

AFGL-TR-01-0069  NL 


O 

01 

o 


< 


Approximate  Analytic  Solutions  for  the  Primary 
Auroral  Electron  Flux  and  Related  Quantities 

JOHN  R.  JASPERSE 
D.  J.  STRICKLAND 


3  March  1981 


DTIC 

f*  AUG  171981^^ 


Approved  for  public  r diltribution  unlimited. 


SPACE  PHYSICS  DIVISION  PROJECT  2310 

AIR  FORCE  GEOPHYSICS  LABORATORY 

HAMSCOM  APB,  MASSACHUSETTS  01731 

AIR  FORCE  SYSTEMS  COMMAND,  USAF 


81  8  17  053 


This  report  has  been  reviewed  by  the  ESD  Information  Office  (OI)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 


This  technical  report  has  been  reviewed  and 
is  approved  for  publication. 


FOR  THE  COMMANDER 


CWief  Scientist 


Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Technical  Information  Center.  All  others  should  apply  to  the 
National  Technical  Information  Service. 


t  /-  1  <j>  t if  ?  J  J'f  j 
/'i.'o  L  -fKt  -  V  '  £ 


Unclassified 


SECURITY  CL  *SS‘P|CATlON  op  This  *»AOE  I'Whwt  0*»»  fenfer«dj 

REPORT  DOCUMENTATION  PAGE  befcm^o^leSg  torm 

t.  REPORT  nCmOER  [T  govt  ACCF$SICNVa  3  PEC?*?' ft-  r"  <* '  4lOG  NUMBER 


1  REPORT  nlmOEH 

AFGL-TR-81-00B9  ^ 


vp-Ajnaqos 


APPROXIMATE  ANALYTIC  SOLUTIONS  FOR 
THE  PKIMARY  AURORAL  ^T.ECTRON  FLUX 
AND 'RBLATED'QUANTITIES  , 

■/<_ _ <£. _ ; 

- _ '* 

John  R./jasperse 
D.  J. /Strickland*  > — 


Sc  tent  i  fir  ..-I  ri  t  ex  - in: 

"T  performing  ORG.  REPORT  NUMBl 

ERP  No.  730  ✓ _ 

8.  CONTRACT  OR  GRANT  NUMBERfsJ 


10  PPOGPAM  ELEMENT^nTOJt 
APE*  A  WO  PK  UHllAniMBC 


9  PERFORMING  ORGANIZATION  NAME  AND  AOORESS  -•-*  .  --  - 

Air  Force  Geophysics  Laboratory  (PHY)  S’ 

Hanscom  Air  Force  Base  /  J  L  J- 

Massachusetts  01731  V  ' 


tt.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  X 

Air  Force  Geophysics  Laboratory  (PHY)  f  /  /  s 
Hanscom  Air  Force  Base  \f 

Massachusetts  01731 


4TMONITORING  AGENCY  NAME  A  ADDRESSflf  dHlete.nl  from  Controlling  O  Hire)  IS.  SECURITY  CLASS,  (of  this  report) 


f  /  J — QJJJJ27  (  i 

\lO  I  23  ICpItO  1  W 


Unclassified 


)  ~1 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (of  the  mbatrmct  entered  in  Block  20,  i l  different  from  Fedort) 


18  SUPPLEMENTARY  NOTES 


Science  Applications,  Inc.,  Vienna,  VA  22180 


19  KEY  WORDS  (Continue  on  revetae  wide  tl  neceaaery  end  identify  by  block  number) 


Linear  transport  theory 
Auroral  electron  flux 
Auroral  energy  deposition  rate 


ZO.  ABSTRACT  (Continue  on  reverae  aide  1/  neceaaary  end  identify  by  block  number) 

In  this  paper  we  show  that  the  linear  transport  equation  may  he  Solved 
exactly  for  the  primary  auroral  electron  flux  in  plane-parallel  geometry  in  the 
forward  scattering  and  average,  discrete,  energy-loss  approximations.  In 
this  approximation  inelastic  scattering  is  taken  into  account  but  elastic  scatter¬ 
ing  drops  out  and  the  solution  is  an  approximation  to  the  flux  in  the  downward 
hemisphere.  Using  the  multiple  scattering  method,  we  obtain  the  solution  as 
a  finite  sum  of  analytic  functions  of  altitude,  energy,  and  pilch  angle  where  — 


,  I™,  1473  COITION  OF  I  NOV  85  IS  OBSOLETE 


Fnclass  i  ficd 

SECURITY  CLASSIFICATION  OF  TmS  PAGE  (*hen  Pete  Entered 


y  o  h  i>  /v  V 


Preface 


The  authors  wish  to  thank  Dr.  A.  L,  Snyder  for  support  and  encouragement 
throughout  the  course  of  the  work  and  Mr.  N.  Grossbard  who  performed  some 
of  the  calculations  for  this  paper. 


:i 


□  D 


Contents 

1.  INTRODUCTION  7 

2.  BASIC  EQUATION  10 

3.  EQUATION  FOR  THE  PRIMARY  ELECTRON  FLUX  10 

4.  FORWARD  SCATTERING  AND  AVERAGE  DISCRETE 

ENERGY-LOSS  A PPROXIAIATIONS  1 1 

5.  MULTIPLE  SCATTERING  METHOD  14 

6.  MULTIPLE  SCATTERING  SOLUTION  15 

7.  HEMISPHERICALLY  AVERAGED  FLUX,  ENERGY 

DEPOSITION  RATE  AND  IONIZATION  RATE  18 

8.  SOLUTIONS  FOR  SPECIFIC  INCIDENT  ELECTRON 

FLUXES  18 

8.  1  Preliminary  Remarks  18 

8.  2  Unidirectional-Monoenergetic  Incident  Flux  19 

8.3  Isotropic-Maxwellian  Incident  Flux  20 

8.4  Isotropic-Monoenergetic  Incident  Flux  21 

9.  RELATIONSHIP  BETWEEN  THE  ENERGY  DEPOSITION 

RATE  AND  CHAPMAN  AND  J  FUNCTK  XS  22 

9.  1  Preliminary  Remarks  22 

9.2  Generalized  Chapman  and  J  Functions  23 

9.3  Energy  Deposition  Rale  for  Unidirectional  Incident 

Fluxes  24 

9.4  Energy  Deposition  Rate  for  Isotropic  Incident  Fluxes  25 

10.  USE  OF  PSEUDOPARTICLES  TO  APPROXIMATE  THE 

SUMS  25 


tmcsDim  h at 


TlLm> 


Contents 


11.  COMPARISONS  28 

11.  1  Preliminary  Remarks  28 

11.2  Comparisons  for  Isotropic-Monoenergetie  Incident 

Fluxes  29 

11.3  Comparisons  for  Isotropic-Maxwellian  Incident  Fluxes  31 

12.  CONCLUSIONS  33 

REFERENCES  37 


Illustrations 


1.  Hemispherically  Averaged  Partial  Fluxes  for  an  Isotropic 

Maxwellian  with  E  =2  KeV,  Q  =  1  erg/cm^  s  and 

N  =  12  0  s  6 

2.  Sequence  of  Partial  Sums  of  Hemispherically  Averaged 

Partial  Fluxes  for  an  Isotropic-Maxwellian  as  in 

Figure  1  27 

3.  Partial  Ionization  Rates  for  an  Isotropic-Maxwellian  with 

Eq  =  2  KeV,  Qg  =  1  erg/cm  s  and  N  =  12  27 

4.  Sequence  of  Partial  Sums  of  Partial  Ionization  Rates  for  an 

Isotropic-Maxwellian  as  in  Figure  3  28 

5.  Energy  Deposition  Rates  from  the  Analytic,  Range,  and  Rees 

Models  for  10,  5,  and  2  KeV  Isotropic-Monoenergetie 

Sources  Each  Containing  1  erg/cm‘  s  30 

G.  Incident  Maxwellian  Energy  Distributions  for  Characteristic 

Energies  of  5,  2,  and  1  KeV  Each  Containing  1  erg/cm^  s  31 


7.  Energy  Deposition  Rates  from  the  Analytic,  Range,  and 

Strickland  Models  for  the  5,  2,  and  1  KeV  Isotropic- 

Maxwellian  Sources  as  in  Figure  6  32 

8.  Differential  Electron  Fluxes  from  the  Analytic  (N  =  12) 

and  Strickland  Models  at  110  km  for  the  2  KeV 

Isotropic-Maxwellian  Source  34 

9.  Hemispherically  Averaged  Fluxes  from  the  Analytic 

(N  -  12)  and  Strickland  Models  at  110  and  1U0  km  for 
the  2  KeV  Isotropic-Maxwellian  Source 


35 


Approximate  Analytic  Solutions  for  the  Primary 
Auroral  Electron  Flux  and  Related  Quantities 


I.  INTRODl  CHON 

The  problem  of  calculating  the  distr  ibution  function  for  auroral  electrons  may 
be  approached  from  two  points  of  view  ,  one  approach  is  to  neglect  collective 
effects  and  treat  short  range  particle -particle  collisions  using  linear  transport 
theory;  the  other  is  to  neglect  short  range  collisions  and  treat  collective  effects 
using  quasilinear  or  nonlinear  plasma  kinetic  theory.  At  high  altitudes  the  latter 
approach  is  often  used  as  the  plasma  is  essentially  collisionless,  whereas  at  low 
altitudes  (below  about  250  km)  linear  transport  theory  is  generally  used  as  the 
plasma  is  neutral  particle  collision  dominated.  The  question  of  how  to  combine 
the  short  range  and  collective  aspects  of  the  problem  is,  and  should  be,  an  active 
area  of  research.  In  this  paper  we  confine  our  attention  to  the  lower  ionosphere 
and  apply  the  methods  of  linear  transport  theory  to  calculate  the  auroral  primary 
electron  flux. 

There  is  an  extensive  literature  on  the  subject  of  electron  transport  in  the 
auroral  ionosphere.  The  various  approaches  may  be  categorized  as  semi- 
empirical,  range  theoretic,  Fokker-Planck,  Monte  Carlo,  and  transport  theoretic. 

K.arly  work  on  transport  properties  of  KeV  auroral  electrons  comes  from 
Chamberlain.  Estimates  were  made  of  ionization  rates  as  a  function  of  altitude 

(Received  for  publication  2  March  1981 

1.  Chamberlain,  J,  W.  (1961)  Physics  of  Aurora  and  Airglow,  Academic  Press,  N.  V. 


for  energetic  KeV  electrons.  The  rates  were  based  on  deposition  r  esults  of 

Spencer-  who  solved  the  electron  transport  equation  using  the  continuous  energy 
9  3  4 

loss  approximation.”’'  Rees  (a  semi-empirical  approach)  applied  an  energy 
dissipation  function  based  on  laboratory  data  by  Grun.  '  The  forms  of  the  func¬ 
tions  obtained  by  Rees  provided  altitude  profiles  of  the  energy  deposition  and 
ionization  rates  for  monoenergetic  and  energy  distributed  sources  with  various 
pitch  angle  dependences. 

Walt  et  al,  were  the  first  to  pr  ovide  altitude,  energy,  anti  pitch  angle  infor¬ 
mation  on  auroral  electron  fluxes  in  the  KeV  range.  To  do  so  they  obtained  a 

numerical  solution  of  a  Fokker-Planck  equation  which  was  originally  used  to  study 

7 

the  properties  of  electrons  trapped  in  the  radiation  belts.  The  method  assumes 

D 

continuous  energy  loss  and  small  angle  scattering.  Banks  et  al,  joined  together 

Walt's  Fokker-Planck  method  and  a  low  energy  approximate  two-stream  transport 

o 

method  by  Banks  and  Nagy,  '  previously  applied  to  photoelectron  transport.  The 
resulting  equation  was  solved  numerically. 

Beger  et  al,  chose  to  examine  auroral  electron  transport  by  applying  Monte 
Carlo  techniques.  *  ’  ^  Backscatter  yields,  backscatter  spectra,  and  altitude 
profiles  of  the  energy  deposition  rate  are  among  the  transport  quantities  that  were 
calculated  by  this  method.  In  the  first  of  the  two  papers  noted  above,  information 
was  also  given  on  the  lateral  spreading  of  KeV,  narrow  electron  beams. 

2.  Spencer,  L.  V.  (1955)  Theory  of  electron  penetration,  Phys.  Rev.  98:1597. 

3.  Spencer,  L.  V.  (1959)  Knergy  dissipation  by  fast  electrons.  National  Bureau 

of  Standards  Monograph  1.  -  -  - 

4.  Rees,  M.  H.  (1983)  Auroral  ionization  and  excitation  of  incident  energetic 

electrons,  Planet.  Space  Sci.  11: 1209. 

5.  Grun,  A.E.  (1957)  Luniineszenz-photometrische  m-  ssungen  der  energie- 

absorption  im  st  rahlungsfeld  von  elect  ronquellen,  emditnensionaler  fall  in 
lu ft,  Z.  Naturforsch. ,  Ser.  A,  12:89. 

(1.  Walt,  M.,  MacDonald,  W.  M. ,  and  Francis,  W.  K.  (1987)  Penetration  of 
auroral  electrons  into  the  atmosphere.  Physics  of  the  Magnetosphere, 

R.  L.  Carovillano,  J.  F.  McClay,  and  H.  R.  Radoski,  eds. ,  D.  Reidel, 
Dordrecht,  Netherlands,  534. 

7.  MacDonald,  W.  M.  and  Walt,  M.  (1981)  Distribution  function  of  magnetically 

confined  electrons  in  a  scattering  atmosphere,  Ann.  Phys.  1  5:44. 

8.  Banks,  P.  M.,  Chappell,  C.R.,  and  Nagy,  A.F.  (1974)  A  new  model  for  the 

interaction  of  auroral  electrons  with  the  atmosphere:  Spectral  degradation 
backscatter,  optical  emission,  and  ionization,  J.  Geophys.  Res.  79:1459. 

9.  Banks,  P.  M.  and  Nagy,  A.F.  (1970)  Concerning  the  influence  of  elastic 

scattering  upon  photoelectron  transport  and  escape,  J.  Geophys.  Res. 
75:1902. 

10.  Berger,  M.J.,  Seltzer,  S.  M.  ,  and  Maeda,  K.  (1970)  Knergy  deposition  by 
auroral  electrons  in  the  atmosphere,  J.  Atmos.  Terr.  Phys.  32:1015. 

11.  Berger,  M.J.,  Seltzer,  S.M.,  and  Maeda,  K.  (1974)  Some  new  results  on 
electron  transport  in  the  atmosphere,  J.  Atmos.  Terr.  Phys.  30:591. 


Ther  o  ar  e  at  least  three  transport  models  utilizing  the  linear  transport 

equation  which  have  rerent ly  been  applied  <>  auroral  studies.  These  are  the 

fa  i  14 

models  of  Strickland  et  al,  “  Mantas  '  and  Stamnes.  Ail  three  models  give  a 
detailed  description  of  elastic  scattering  and  ullo.v  for  discrete  energy  loss.  Dif¬ 
ferences  arise  >n  the  representation  of  the  flux  within  the  collision  integral  and  in 

the  method  of  integration  over  depth.  Strickland  et  al,  allow  the  flux  to  vary 

19 

quadrat  ically  in  fn  K  and  .i  nearly  in  u  within  ;tny  given  T,  w  cell.  "  The  quadratic 
dependence  was  introduced  because  of  energy  conservation  problems  for  a  linear 
dependence  when  treating  energetic  fluxes  above  several  KeV.  The  integration 
over  depth  is  carried  ou‘  by  either  a  finite  difference  method  (second  order 
pred  etoi -corrector)  or  an  eigenvalue  method.  The  latter  approach  was  found  to 
be  mu'h  faster  and  more  accurate.  Mantas  allows  the  flux  to  vary  linearly  in 
both  T  and  u  within  any  K,  u  cell.  ^  He  also  uses  a  linear  dependence  within  a 
given  '  cell  which  leads  to  the  standard  finite  difference  expression  for  the  first 
order  z  derivatives.  Stamnes  treats  the  y  dependence  of  the  problem  by  the  dis¬ 
crete  ordina'e  method  and  considers  the  flux  to  he  constant  within  a  given  If  cell 

14 

(commonly  called  the  multigroup  approximation).  Tor  a  good  discussion  of  the 

la  1 0 

discrete  ordinate  and  multigroup  methods  see  Davison.  ’  Like  Strickland  et  al, 

14 

Stamnes  uses  an  eigenvalue  technique  to  carry  out  the  integration  over  depth. 

The  discrete  ordinate  method  for  treating  the  angular  dependence  of  the  electron 
flux  converges  poorly  for  highly  anisotropic  scattering  kernels.  Therefore,  the 
angular  dependence  of  the  results  of  Stamnes  at  high  energies  in  the  backward 
direction  are  questionable.  All  of  the  transport  methods  lead  to  a  truncated 
matrix  equation  which  is  solved  numerically. 

7 ne  new  feature  of  the  work  presented  by  us  in  this  paper  is  that  we  give 
approximate  analytic  solutions  for  the  primary  auroral  electron  flux  in  the  down¬ 
ward  hemisphere  and  related  quantities.  All  previous  work  on  this  problem  has 
resulted  in  numerical  solutions.  The  analytic  results,  though  approximate,  are 
useful  in  giving  insight  into  the  physics  of  auroral  electron  precipitation  and  in 
carrying  out  further  analytical  studies  of  auroral  phenomena.  For  example,  in 
work  currently  in  progress  using  the  multiple  scattering  method,  we  show  that 
elastic  scattering  can  be  included  rigorously  and  that  the  solution  given  in  this 

12.  Strickland,  D.J.,  Book,  D.  L. ,  Coffey,  T.P.,  and  Tedder,  J.A.  (197(1) 

Transport  equation  techniques  for  the  deposition  of  auroral  electrons, 

J.  Geophys.  lies.  81:2755. 

13.  Mantas,  G.  P.  (1975)  Theory  of  photoelectron  thermalization  and  transport 

in  the  ionosphere.  Planet.  Spare  Sci.  23:337. 

14.  Stamnes,  K.  (1978)  A  theoretical  investigation  of  the  interaction  of  auroral 

electrons  with  the  atmosphere,  Ph.  D.  thesis,  I'niversity  of  Colorado, 
Colorado. 

15.  Davison,  B.  (1957)  Neutron  Transport  Theory,  Oxford  Press,  London. 
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paper-  is  just  the  leading  teem  in  an  iterative  sequence  of  solutions  which  treats 
the  complete  problem  including  secondary  electron  production.  In  work  also  in 
progress  the  method  presented  in  this  paper-  is  used  to  solve  the  coupled  proton- 
hydrogen  precipitation  problem  for  the  proton  aurora.  These  Solutions  for  the 
proton  and  hydrogen  flux  and  related  quantities  are  similar  to  the  ones  presented 
here  for  the  electron  precipitation  problem. 


2.  BASIC  KQl.  ATION 


For  a  one  constituent  atmosphere  the  linear  transport  equation  for  energetic 
electrons  in  plane-parallel  geometry  is 


—  ■  n(z)Q((K) 


<t>  ( / 


,  K,  nlzlQ^K)  J  cl K '  dy'  |{(K'u'  -  Kyi  $(/,  IT,  y ') 


( 1 ) 


where  <f>  is  the  electron  flux,  n  is  the  neutral  density,  and 

K(KV  -  £  adl-iV  .  F.y)  .  (21 

1  j 

Here,  /  is  the  altitude,  K  the  particle  energy,  y  the  cosine  of  the  angle  between 
the  particle  velocity  and  the  positive  z-axis,  and  o  the  differential  cross  section 
per  unit  energy  range  (referred  to  hereafter  as  the  differential  cross  section)  for 
elect ron-neut ral  particle  scattering  where  the  neutral  particle  makes  a  transition 
from  the  ground  state  to  the  final  state  i.  The  total  cross  section,  Q  (Hi,  is  re¬ 
lated  to  the  differential  cross  section  by  the  formula 

Q,(K)  =  2*  J  di-T  dy'  Y,  ‘^Cly  -  MV)  .  til 

t.  K(,M  ATION  FOR  TIIK  I’RIMAIO  F.l.KCTHON  KIT  \ 

The  differential  cross  sections  describe  elastic,  excitation,  and  tom/ation- 
type  collisions  between  the  electrons  and  the  neutral  particles.  In  ioni  rat  toti-tvpe 
collisions  additional  (secondary)  electrons  are  produced.  The  electron  flux  in 
Kq.  (1)  includes  both  primary  and  secondary  electrons.  The  secondary  election 
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sport  rum  is  highly  peaked  at  low  energies  with  very  few  secondary  elec- 

t  funs  at  higher  energies  ('>K|nln>-  1,1  In  the  auroral  ionosphere  h'mm  is  somewhat 
arbitrary  but  is  on  the  order  of  1  2  KeV  for  primary  electrons  ranging  from 
1  to  (i 0  KeV.  If  we  restrict  the  energy  range  of  Kq.  (1)  to  high  energies,  we  may 
approximate  the  electron  flux  by  neglecting  the  secondaries  eompaied  to  the 
primaries;  that  is, 


4>  .  fo 

primaries 


K  ^  i: 


(4) 


and  use  a  differential  i*roSS  section  for  ioiii /at  ion-type  collisions  that  describe:- 
only  the  scattering  of  the  primary  electrons.  Viewing  the  ioni/ation-type  colli¬ 
sions  in  this  wav  we  may  approximate  the  differential  et'oss  sections  for  KeV 
auroral  electrons  bv 


>  (KV  -  Ku>  (3 " 1  i :■> *( i -  i:>  ^.(K'.m'.u)  , 


c>> 


,Vxku:v  .  Koi  '  a,^'1yi>xk(i:')6(i:1  -  (i:  -  wk)l  iYxk(u'.K'  .  d 


-  K u)  '  (2rr)'1n|l,)  <i\  O  p J'/V.u)  . 


lliMi',  ihf  subsoripts  t',  in,  nnti  i  i  ■  i '  ft '  i  to  olast  ir,  oxoitation,  anil  toni /at  ii  in -i  ypo 
.  ollisions,  i  i'spoi'tivi'lv.  6  is  tin-  Dii  ao  il.-lta  funotion,  tin-  !'  funotiotis  aro  sym- 
niririi'  mu’  atni  a  anil  air  normali/od  to  otio,  thi'  total  oross  sort  ion  for  oaoli 
proross  is  ilonotmi  In'  Q,  anti  tho  supo  rso  n  pt  p  rofors  to  primary  olfi-t  rolls.  \\  .. 


•I, 


■lnsti. 


oxoitnt  ion.  am!  r  -ti i 


1 


is  tho  oxoitation  onorgy.  l  or  a  ilisrussi 
t  ion  tvpo  rross  sooiions  soo  St  riokland  o 

Tilt'  oquat  ion  fo  r  tho  p  rim  a  ry  Hoot  poti  flux  is  I  iq.  (1)  with  'In'  ii  i  ffo  rotit  i  a  1 
oross  soot  ions  givon  by  Kqs.  (a),  ('>).  and  (7),  and  is  valid  for  1.  >  1. 


i  K)R\OKi>  scvitfrim;  \m»  wkh  \<;k  kisckktk  km  k(;v 

I  .OSS  \ITRO\IM  VTIONS 

In  tho  auroral  rogioti  two  faots  about  tho  difforontial  oross  sootiotis  at  high 
onorgios  aro  appa  rotit.  Khov  aro  that  1)  tho  tii  ffo  rout  i  a  1  oross  soil  ions  an  lughl\ 
poakt'd  in  tho  forward  dirootioii,  iho  oxoitaiioii.  and  ionization  tvpo  ooUisioii.-'  bring 

1U.  ( ipal,  (  .Ik.  I’otorsoti,  \V.  K.  .  and  Hoaty.  li.C  .  <l!‘7l>  Moasuroitiotits  of 

soooiida  rv-oli'ot  roll  spool  ra  proiiuood  hv  olool  roil  mipai't  tom/at  i- >; .  ‘  a 
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r 


even  more  highly  peaked  than  the  elastic  collisions,  and  2)  the  average  energy 
loss  per  inelastic  collision  (excitation  and  ionization)  is  a  weak  function  of  the  in 
cident  electron  energy.  This  suggests  that  for  inelastic  collisions  we  make  the 
forward-scattering  and  average  discrete  energy-loss  approximations: 

I  °exk<EV 
k 


f 

where  W  is  the  average  primary  electron  energy  lost  per  inelastic  collision 
defined  by 


Z  Qexk<K,)  )  6lE'  "  (E  +  w>!  -  «> 

k 


(8) 


£  Q^fE')]  6[E’  -  (E  +  W)]  6(p'  -  p)  ,  (9) 
t 


W  = 


[IWk*exk(E>  +  Z/dE,(E-E’>0!,P)(E 


(E,  E') 


-1 


[I  Qexk<E>  +  l/dE'oJP>(E,E.)' 

L  k  t 


(10) 


W  as  a  function  of  E  for  the  auroral  cross  sections  used  by  Strickland  et  al,  is 

32,  30,  42  and  48  in  eV  for  1,  3,  10  and  30  KeV  primary  electron  energies, 

12 

respectively.  In  the  average  discrete  energy-loss  approximation  \Y  is  assumed 
constant.  For  elastic  cross  sections  a  screened  Rutherford  cross  section  is 
often  used: 


ap<EV  -  Ep)  -  (2tr)‘1Qe(E')6(EI  -  E)  p (n.p'.p)  . 


(11) 


p(n, p) 


2p(l  +  r))(l  +  2q  -  p  'p ) 


t(l  +  2n  -  p'p)2  ~  (l  -  p,2H1  -  u2)\ 3/2 


(12) 


where  i)  is  a  function  of  energy.  For  graphs  of  p  as  a  funtion  of  q,  p1  and  p,  and 

12 

n  as  a  function  of  E  see  Strickland  et  al. 


Inserting  Kq.  (ii)  through  Kq.  (12)  into  Kq.  (1)  and  transforming  Kq.  (1)  ti 
an  equation  m  terms  of  the  optical  depth,  T,  where 

dr  -n(z)  Q(K)  dz  ,  ( 


we  obtain 


4 >(t,  K,  u ) 


Qd:  -r  \\  ) 


<Mt (z,  i;  *-  w),  i;  •  w. 


Q  <K) 


<mt,  i:, q )  - 


J  'la  '  p(q,  a a^  4><T,  II,  a  'I 


where  Q  and  Q(i  denote  the  total  inelastic-  and  total  elastic  cross  sections,  respec 
tively,  and  t  denotes  t(z,  K).  The  ratio  Q(|;  •  W)  'Q(E)  is  near  1  and  the  ratio 

QJKl  Q(K),  as  given  by  Strickland  et  a),  is  0.84,  0.82,  0,45  and  0.27  for  1,  5‘, 

*'  1  ■> 

10  and  .20  KeV  primary  electron  energies,  respectively.  ~  (in  the  right-hand  sid 

of  Kq.  (14)  we  may  view  the  second  and  third  terms  as  producing  a  correction  to 
the  first  term.  This  correction  is  not  necessarily  small,  l  or  example,  we  know 
that  in  the  backward  hemisphere  elastic  scattering  is  the  dominate  process  and 
contributes  the  most  to  the  backseattered  flux.  The  solution  here  is  very  sensi¬ 
tive  to  the  ratio  Q  and  to  the  precise  shape  of  p(ij,  p ’,  (1 ).  However,  in  the 
downward  hemisphere  we  expect  the  first  term  on  the  right-hand  side  of  Kq.  (14) 
to  produce  a  contribution  to  <1- which  is  comparable  to  or  larger  than  the  contribu¬ 
tion  made  by  the  second  and  third  terms.  In  this  paper  we  seek  an  approximation 
to  the  flux  in  the  downward  hemisphere.  An  equation  for  such  an  approximate 
solution  is  obtained  from  Kq.  (14)  by  taking  the  limit  of  the  function  p  as  K  -  x 
(t?  -  0).  It  can  be  shown  that 


pin,  U  q  >  —  fliu  1  -  q) 

17-0 


In  the  high-energy  approximation  elastic  scattering  becomes  highly  peaked  in 
the  forward  direction  and  drops  out  of  Kq.  (14).  Wo  obtain 


;  %.  -  a 


*(T(Z,  E),  !•;,„)  <£(T(z,  K  4  W),  K  +  W,q)  .  (lli) 


Equation  <  1  (>)  is  a  partial  differential-difference  equation  for  the  primary  election 

flux  now  denoted  by  <$,  as  a  function  of  r,  K  and  p,  approximately  valid  for 

0  ■<  T  <  oe.  E  «  K  <  op.  and  -1  «?  u  *-  0. 
nun 


5.  Ml  LTIP1.K  SCATTERING  MKTIIOI) 

The  multiple  scattering  method  may  be  applied  in  principle  to  any  linear 
transport  equation  satisfying  certain  general  requirements.  For  a  discussion  of 

the  multiple  Scattering  method  applied  to  elastic  scattering  see  Goudsmit  and 

17  18  19 

Saunderson  and  Wang  and  Guth,  to  inelastic  scattering  see  I-'ano,  and  to 

20 

both  types  of  scattering  see  Case  and  Zweifel.  “  In  our  case  w  e  w  ish  to  apply  it 
to  Kq.  (lti)  which  contains  only  inelastic  scattering.  To  do  so  we  rewrite 
Kq.  (lti)  as 


M  §7  '  1 )  *  & 


(17) 


where  oC  is  the  linear  operator  that  shifts  the  energy  variable  of  by  W  and 
multiplies  by  Q(K  +  \V)/Q(K).  The  multiple  scattering  method  consists  of  writing 


00 

*(r,  E,p)  -  Y,  ^n(T-  •  08) 

n  0 

where  the  <f> ^  satisfy  the  infinite  system  of  equations 


0 


(19) 


(20) 


17.  Goudsmit,  S.  and  Saunderson,  J.  L.  ( 1940)  Multiple  scattering  of  elect rons, 

Phys,  Hev.  57:24. 

18.  Wang,  M.C.  and  Guth,  E.  (1951)  On  the  theory  of  multiple  scattering,  par¬ 

ticularly  of  changed  particles,  Phys.  Rev.  84:1092. 

19.  Fano,  U.  (1953)  Degradation  and  range  straggling  of  high-energy  radiations, 

Phys,  Rev.  92:330. 

20.  Case,  K.  M.  and  Zweifel,  P,  L.  (1967)  Linear  Transport  Theory,  Addison- 

Wesley,  Reading,  48. 
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and  the  solution  6  in  terms  of  <4  ,  for  n  1  is 

vn  ^n-l 


^n(T,  K,  n)  =  -  exp  (t/jj)  J  dtp  1  b1 1  <J>n_1(b11t,  K  +  W,  u)  exp  (-l/u)  ,  (2t>) 

o 

for  -1  s  p  <0.  For  0  <:  n  ■<  1  the  are  zero.  The  complete  solution  is 

00 

$(t,  K,  m)  =  <5>(F,  p)  exp  (T/ji)  +  £  $(10  +  nW,  p)  b  (K) 

n- 1 


x  [IIn(T/w,  10)  -  iyo,  10)  exp  (t /jj ) ]  , 

(27) 

for  -1  <  ji 

tions  is 

<.  0,  and  $  -  0  for  0?u  <  1.  The  recursion  relation  for  the  11 

funr  - 

exp  (- 

x)  Hn+1(x,  10)  -  Hn+1(0.  K)  - 

X 

-  J  dx'  exp  (-x')tHn(bux»,  E  +  W)  -  11^(0,  K  +  \V)  exp  (bnx')l 

O 

.  (28) 

An  explicit  formula  for  the  general  H  function  has  been  worked  out  by  mathemat¬ 
ical  induction  but,  as  it  is  quite  lengthy,  only  the  first  three  functions  w  ill  be 
given  here: 

J!t<x,  1 

10)  a  j  j  exp  (bj  jX) 

(20) 

1 

K)  a22a21  oxp  (b21x)  ‘  :I22:| 1 1  oxp  <bux>  . 

(20) 

H.<(x,  1 

K)  =  a:i3a32a3l  oxp  (b:ux)  '  a:i.rn:i2aii  oxp  (bux> 

‘  a33a22  '21  °xp  (b21x)  +  a33a22a  1 1  °Xp  (bllx)  ' 

(21) 

In  the  above  formulae  we  have  introduced  the  following  definitions: 

1) 

nm 

bnni(K)  -  Q0-:  +  n\V)  '(■)(  10  +  (m  -  1)\V)  , 

(2  2) 

Hi 


(33) 


a  =  a  (E)  =  (1  -  b  <E)]  _1  . 

nm  nm  nm 

Equation  (27)  may  be  rewritten  in  an  alternative  form  as 

oo  n 

$(t,  E,p)  =  4>(E,/j)  exp  (t/h)  +  £  $(E  +  nW,p)  £  A 

n=  1  1  =  1 


X  (exp  (Bj  t/^)  -  exp  (t/^))  , 


(34) 


where 

Bf(E)  =  Bi  =  bM(E)  ,  (3  5) 

and  the  first  few  Anf  =  An^(E)  are 


A11 

=  +bUau  . 

(3fi) 

A22 

=  +  b21a22a21  ’ 

(37) 

A21 

=  "  b2  la22a  1 1  - 

(38) 

A33 

=  +  b31a33a32a31  ' 

(39) 

A32 

b31a33a22a21  ’ 

(40) 

A31 

“  b31^+  a33a32all  "  a33a22all^ 

(41) 

Several  comments  about  this  solution  are  in  order.  The  primary  electron 
flux,  $,  is  an  exact  solution  to  Eq.  (10)  which  is  in  turn  an  approximation  to 
Eq.  (1).  The  fact  that  it  solves  Eq.  (1(5)  exactly  can  be  verified  by  direct  sub¬ 
stitution.  The  value  of  $  depends  linearly  on  its  boundary  value;  doubling  the 
incident  flux  doubles  the  flux  at  all  altitudes.  The  solution  yields  no  backscattered 
primary  electron  flux.  The  reason  for  this  is  that  in  the  forward  scattering  ap¬ 
proximation  elastic  scattering  drops  out  and  no  particles  are  scattered  into  the 
backward  direction. 
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7.  HEMISPHERICALLY  AVERAGED  FLUX,  ENERGY  DEPOSITION 
RATE  AND  IONIZATION  RATE 

The  hemispherical  electron  flux,  is  defined  as  the  average  value  of  the 

electron  flux  over  the  downward  hemisphere  and  is 


$h(t,  E) 


o 

2ir  J  dp  $(t,  E,  p)/2tr  .  » 

-1 


(42) 


The  energy  deposition  rate,  )?  is  defined  as  the  energy  deposited  per  unit  vol¬ 
ume  per  unit  time  by  the  precipitating  electrons,  and  is  given  in  plane-parallel 
geometry  by 


oo  +1 

rjj,(z)  =  2»  J  dE  J  dpuEn(z)Q(E)  4>(t,  E,g)  .  (43) 

o  -1 

The  ionization  rate,  rj . ,  is  defined  as  the  number  of  electron-ion  pairs  produced 

per  unit  volume  per  unit  time.  This  can  be  found  to  a  good  approximation  from 

2 1 

nF  by  using  Bethe's  formula 

n.(z)  nKU)/Ks  .  <44> 

where  E  is  a  constant.  A  good  value  for  E  for  auroral  electrons  has  been 
s  ”  s 

found  to  be  34  eV.  In  this  way  we  find  an  approximate  formula  for  rj;  without 
solving  explicitly  for  the  secondary  electron  flux. 

8.  SOU  TIONS  FOR  SPECIFIC.  INCIDENT  ELECTRON  FLUXES 
8.1  Preliminary  Remark* 

In  this  section  we  give  analytic  expressions  for  the  quantities  defined  in 

Section  7  for  several  specific  Forms  of  the  incident  electron  flux.  Before  we  do 

so,  we  remind  the  reader  that  our  analytic  solution  for  the  electron  flux  is  only 

valid  for  E  >  E  .  .  However,  in  the  energy  integral  from  0  to  oc  which  defines 
m  tn  ’  ”•  " 

we  note  that  the  contribution  to  the  integral  below  E  .  only  amounts  to  a  few 
E’  °  min 

percent  for  auroral  electrons. 

21.  Bethe,  H.A.  (1933)  Iianbueh  der  Physik,  Verlag  Julius  Springer,  Berlin 
24:491. 
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8.2  l)nidirectional-Moiioener®etic  Incident  Flux 

For  a  unidirectional-monoenorgotic  electron  flux  incident  at  the  top  of  the 
ionosphere  the  boundary  condition  is 


<I>  (E,  p)  = 


6<E  -  EQ)6(p  +po) 


(45) 


where  Qg  is  the  total  energy  flux  in  the  downward  direction,  E()  is  the  electron 
energy,  and  po  is  the  cosine  of  the  angle  of  incidence,  x,  where  0  ■s  \  <  *'2. 
The  electron  flux  is 


4>(T,  E,p)  = 


6(p 


“o' 


{' 


6(E  -  Eo)  exp 


+ 


Y  6(E  +  nW  -  E  )b  .  (E)  til  (T/p,  E)  -  II  (0,  E)  exp  (T/p)] 
*•—>  o  n  i  n  n 

n=l 


The  hemispherical  flux  is 


<fH<T,  K)  = 


6(E  -  E  )  exp  (-T/p  ) 

o  o 


(4ti) 


+  V  6(E  -  (E  -  nW))b  ,  [H  (-r/p  ,  E)  -  II  (0,  E)  exp  (-T/p  )]) 
o  nj.  n  o  n  o  I 

n  =  1 


The  energy  deposition  rate  is 


*E<Z>  =  \~)  Q(Vn(z) 

■  O  # 


|exp  (-T(z,  EQ)/po>  +  y 
^  n  =  1 


(1  -  nW/E  )  b  ,  ( E  -  nW) 
o  n  I  o 


X  [In(-T(z,  Eo  -  nW)/pn,  E^  -  nW)  -  HJO,  En  -  nW) 


X  exp  ( -t(z,  Eq  -  nW)/po))j> 


(48) 
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where  the  I  functions  are 


In(x,  E)  = 


3 

3x 


Hn(x,  E) 


(49) 


It  should  be  pointed  out  here  that  the  above  sums  do  not  extend  to  ao  but  to  a  max¬ 
imum  value,  N,  given  by  NW  =  K  -  K  .  . 

’  B  -  o  min 


8.3  Isotropic- Maxwellian  Incident  Flux 

For  an  iSotropic-Maxwellian  electron  flux  incident  at  the  top  of  the  ionosphere 
the  boundary  condition  is 


$(E, p) 


E  exp  (-E/E^) 


for  -l<(i<  0 


(50) 


where  Q  is  the  total  energy  flux  in  the  downward  direction  and  E  is  the  charac- 
s  ”•  o 

ter  istic  energy.  The  resulting  electron  flux  is  given  by 


<5>(r,  E,p) 


E  exp  <-E/E'o)  exp  (T/p)  + 


x 

n  =  1 


(E  +  nW) 


X  exp  (-  (E  +  nW)/Eo)  b 


nl  [Hn(x /p,  E)  -  Hn(0,  E)  exp  (r/M)J 


(51) 


with  $  =  0  for  0  ^  p  -  1,  The  hemispherical  flux  is 

00 

E  exp  (-E/Eo)  E2(t)  +  £  (E  +  nW) 
n=  1 

X  exp  (-(E  +  nW)/E  )  b  .  [K  (t,  E)  -  H  (0,  E)  E„(r)l 
o  nl  n  n  l 

where  Egfx)  is  the  second  exponential  integral  function  defined  by 

CO 

E2<x)  /  dt  t  ^  exp  (-xt) 

1 


(52) 


(53) 


20 


iiikr  :***•-  rih  daniiti 


■y- 


The  first  three  K-functions  are 


K^x,  E)  a11E;,(b11x) 


(54) 


K9(x,  E)  =  a92a2jE2(b21x)  -  a.>2an  E2(t>2  lx)  ’ 


(55) 


K3(x,  E)  -  a33a32a3iE2(b31x^  '  a33a32allE2(bl lx) 

'  a33a22a21K2<b21x)  +  a33a22a 1 1  L2(bl lx*  * 


(5(i) 


The  general  rule  for  the  K  functions  is  the  same  as  that  for  the  H  functions  with 
the  exponential  function  replaced  by  the  second  exponential  integral  function. 
The  energy  deposition  rate  is 


Q. 

rlE(z)  =  |  ,,3 

sl 

o  i 


00 

n(z)  J  dEEQ(E)  |e  exp  (-E/Eo)  E2(t) 


E 

min 


+  Y.  +  nW)  exp  <-(E  +  nW)/Eo)bnl[Ln(T,  E)  -  Hn(0,  E)E2(t)) 
n=  1 


(57) 


The  rule  for  the  L  functions  is  the  same  as  that  for  the  I  functions  except  that  the 

exponential  function  is  replaced  by  the  second  exponential  integral  function.  In 

performing  the  integral  over  the  total  inelastic  cross  section,  Q,  the  integral 

may  be  truncated  at  some  E  .  When  this  is  done  the  infinite  sum  terminates 

max 

at  a  maximum  value,  N,  given  by  NW  -  E  -  E  .  . 

’  '  b  J  max  nun 


8.4  laotropic-Monoenergetii  Incident  Flux 

For  an  incident  isotropic -monoenergetic  flux  the  boundary  condition  is 


<t(E,p)  = 


for 


-1  £  p  <  0 


(58) 


where  Q  is  the  total  energy  flux  in  the  downward  direction  and  E  is  the  electron 
®  o 

energy.  The  expression  for  4>  is 
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$(T,  E ,n)  =  |fi(E  -  Eo)  exp  (t/u)  +  Y  ME  +  nW  -  Eq) 


X  [Hn(T/^,  E)  -  Hn(0.  E)  exp  <r/M)]j 


with  4>  0  for  0-?  u  s  1. 


9.  RELATIONSHIP  RETWEEN  THE  ENERGY  DEPOSITION  RATE 
AM)  CH  APM  AN  AND  J  El  NOTIONS 


9.1  Pri'limiiiar)  Remarks 

Chapman  functions  are  functions  of  altitude  which  arise  in  the  Chapman  theory 
of  the  daytime  ionosphere.  The  functions  gives  the  ionization  rate,  which  is  re¬ 
lated  to  the  energy  deposition  rate,  for  a  single  Chapman  layer.  The  ionization 
rate  profile  is  produced  by  (he  photoionization  of  the  neutral  particle  gas  by 
ionizing  electromagnetic  radiation  from  the  sun.  In  the  auroral  zone  at  night, 
the  ionized  layer  is  produced  by  a  different  process;  the  precipitation  of  energetic 
electrons  incident  at  the  top  of  the  atmosphere.  It  has  been  known  for  some  time 
that  the  auroral  ionization  rate  profile  (or  energy  deposition  rate  profile)  is  much 
more  highly  peaked  as  a  function  of  altitude  than  is  a  single  Chapman  function. 

In  this  section  we  establish  the  connection  between  the  Chapman  functions  and 
the  auroral  ionization  rate,  specifically  we  show  that  the  solution  for  the  energy 
deposition  rate  in  the  auroral  ionosphere  is  a  superposition  of  Chapman  functions 
for  unidirectional  incident  electron  fluxes,  and  a  superposition  of  what  we  call 
J  functions  for  isotropic  incident  electron  fluxes.  In  order  to  do  this  wo  use  the 
alternative  expression  for  r)j.  we  obtain  when  Eq.  (34)  is  substituted  mte  Eq.  (43). 


r)j.(z)  2*  n(z)  J  dE  J  du  EQ(1.)  j<t(E,p)  exp  (T(z,  E)  w) 

E  .  - 1  ^ 

nun 


£  $(E  +  n\V,M)  Y  Ant{K) 


[B  (K)  exp  (B  (i;)T(z,  -  exp  (t(z,  K)  w) |  ( 


9.2  Generalized  Chapman  and  J  Function* 


The  Chapman  functions  as  introduced  by  Chapman  is 
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Ch(z)  =  n(z)  exp 


(Gl) 


where  f  is  a  constant  and  is  the  cosine  of  the  angle  of  incidence,  where 
0  s  x  <  't/s.  Let  us  now  define  a  generalized  Chapman  function  of  three  argu¬ 
ments 


Ch  2,  —  f(K)  =  n(z)  exp 
1  P  ' 


/  dy  n(y) 


(02) 


where  Os  z  <«,  -1  <  ji  <  0,  and  f  is  a  positive  function  of  E,  where  0  <  E  <  ®. 
By  analogy  with  the  Chapman  function  let  us  introduce  a  J  function 


J(z)  s  J  dy  n(z)  exp 


•j  /  dy  n(v) 


(63) 


and  a  generalized  J  function 


J(z,  f(E))  s  J  d(j  n(z)  exp 


-1 


^  /  dv  n(y) 
P  J 


(04) 


where  Os  z  <  oo,  and  f  is  a  positive  function  of  E,  where  Os  E  <  oc.  Transform¬ 


ing  Eq.  (04)  to  an  integral  on  t,  where  t  =  -p  *,  we  see  that 


J(z,  f(E))  =  n(z)  E2 


f(E)  J  dy  n(y) 

z 


(05) 


where  E9  is  defined  by  Eq.  (53). 


22.  Chapman,  S.  (1931)  The  absorption  and  dissociative  or  ionizing  effect  of 
monochromatic  radiation  in  an  atmosphere  on  a  rotating  Earth,  Pipe. 
Phys.  Soc.  43:20. 
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9.3  Energy  Deposition  Hate  for  l  nidirectional  Incident  Fluxes 

For  a  unidirectional  incident  flux  (p  cos  x.,  Us  X  <  n/2)  the  energy  depo¬ 
sition  rate  given  by  F.q .  ((>0)  is 

E 

max 

nF( z)  --  2»  J  dEKQ(K) 

E 

min 


N  n 

+  I  I  *(E  +  n\V)An#(K) 
n- 1  /  1 


. 

\  / 

Bf(K)Ch  L 

B((K)Q(E )j  -  Ch  f 

(tin) 

Thus,  the  exact  solution  of  Eq.  (1(>)  yields  the  energy  deposition  rate  as  an  energy 
integral  over  a  finite  superposition  of  generalized  Chapman  functions. 

The  explicit  solution  for  a  unidi reel ional-monoeno igetic  incident  flux  - 


nF(z) 


r-. 


N  n 

+  X  I  <!  -  n\V/En)B((Eo  -  nWlA^  -  n\V) 
n  - 1  f  - 1 
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VK, 


-  n\V)  Ch  7. 


-  — ,  B,  (FI  -  n\V )  Q( E  -  n\V)  j 
a  f  o  o  / 
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This  explic  it  solution  reduces  to  a  finite  superposition  of  generalized  Chapman 
functions  with  shifted  energy  arguments. 


9.4  Knrr»>  l)«*|H>sili<>n  Kali-  for  Isotropic  Incident  Hiiws 

For  :in  isotropic  incident  flux  the  energy  deposition  rate  is 


N  n 


nK(/)  i>tr  J  di:  i:q(k) /<t(K)J(z,gu:))  •  £  £  *<i-:  •*  n\v)  An(  (K) 


n  1  f  1 


V  [Bf  (K)J(z,  l?f(i:)Q(K))  -  J(/,Q(i;))] 


(lib) 


Here,  the  exact  solution  of  Fq.  (Ill)  yields  an  energy  deposition  rate  as  an  energy 
integral  over  a  finite  sum  of  generalized  J  functions. 

For  an  isot  ropic-Monoonerget  ie  incident  flux  the  enery  integral  can  be  done 
explicitly  and  for  an  isot  ropic-Maxwellian  incident  flux  the  energy  integral  can  be 
done  by  quadrature. 


10.  I  SK  OK  PS K l  1)01’  VRTICl.KS  TO  VITROMM  \TK.  TIIK  SI  MS 

In  the  formulae  of  Sections  0  through  0.  W  is  a  number  on  the  order  of  40  eV 
and  the  primary  auroral  electron  energies  range  from  about  i  to  00  KeV.  As  a 
result,  the  sums  in  the  formulas  extend  to  as  many  as  1500  terms  representing  as 
many  as  1500  scatterings  before  an  energetic  electron  loses  all  its  energy.  How¬ 
ever,  it  turns  out  that  a  good  approximation  to  these  large  sums  may  be  obtained 
by  introducing  the  notion  of  pseudopartieles  which,  in  turn,  allows  us  to  take 
many  fewer  terms  in  each  sum.  A  pseudoparticle  is  a  particle  which  as  a  cross 
section  \V  -'W  times  smaller  than  the  real  particle  but  has  an  average  energy  loss 
per  inelastic  collision  \V  W  times  greater,  such  that  the  product  of  the  two  re¬ 
mains  the  same.  With  ;tiis  approximation  a  good  answer  is  obtained  for  the  sums 
given  in  Sections  t>  through  P  with  many  fewer  terms  in  each  sum. 

In  order  to  introduce  pseudopartieles  into  the  equations  of  Sections  il  through 
P  we  simply  replace  Q  by  and  W  by  \V^,  that  is, 

Q(H)  —  Q  (K)  =  (W  W  )  Q(E)  ,  (IIP) 

P  P 

W  -  W  —  (W  /\V)W  .  (70) 

P  P 


We  have  examined  the  convergence  of  the  pseudopart iele  method  as  we 
approximate  the  sums  by  a  smaller  and  smaller  number,  N,  of  pseudoseatterings 


spanning  the  energy  range  from  ^  to  K  .  We  have  found  that  a  good  answer 

can  be  obtained  with  a  surpr  isingly  small  number  of  pseudoseatter mgs.  For 
example,  for  an  tsot  ropie-Maxw  ellian  incident  flux  with  H  ranging  between  1  anti 
10  KeV  only  10  to  15  pseudoseatterings  spanning  the  energy  range  are  needeti  to 
obtain  answers  to  within  a  few  percent  of  the  full  Solution  it  most  altitudes  of 
interest. 

In  Figure  1  we  show  the  hemispherieally  averaged  partial  fluxes,  Ojjn,  for  a 

sequence  of  n  values  for  an  isot  ropie-Maxw  ellian  incident  flux  w  ith  K  2  KeV 

2  11 
and  =  1  erg,  cm“s  for  the  ease  where  N  12.  That  is  to  say,  12  pseudoseat- 

terings  spanning  the  energy  range  of  interest  from  K  -  11  Ke\’  to 

K  1  KeV.  In  Figure  2  we  illustrate  the  rate  of  convergence  of  <£,,  bv  show  - 

ing  the  sequence  of  partial  sums  of  for  the  same  ease.  In  Figure  J  we  show 

the  partial  ionization  rates  (n  -  n,.  /34)  for  the  same  ease.  Note  that  for  n 
'  in  I.n  to 

up  obtain  a  J  function  profile  and  for  17  (n  :v  1)  we  obtain  functions  with  a  single 
node,  each  of  which  subtracts  from  r).  at  high  altitudes  and  adds  to  rt.  at  low 
altitudes.  In  Figure  4  we  show  the  rate  of  convergence  of  by  show  ing  the 
sequence  of  partial  sums  of  »).  for  the  above  case. 


(Z,E  =  2  5KeV)  (el/crnZ  s  eVsrl 


Figure  1.  Hemispherieally  Averaged  Partial  Fluxes  for 
an  Isotropic  Maxwellian  with  E  =2  KeV.  0=1  erg / 
rm  s  and  N  =  12  °  s 


Figure  2.  Sequence  of  Partial  Sums  of  Ileniispherieallv 
Averaged  Partial  Fluxes  for  an  Isot ropie -Maxwellian  as 
in  Figure  1 
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Figure  3.  Partial  Ionization  Kates  for  an  Isotropie- 
Maxwellian  with  E  =  2  KeV,  Q  =  1  erg/cnms  and 
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Figure  4.  Sequence  of  Partial  Sums  of  Partial  Ioniza¬ 
tion  Hates  for  an  Isot ropic-jllaxwel) ian  as  in  Figure  3 

As  a  general  rule  we  find  that  the  lower  the  altitude  the  more  pseudoscatter- 
ings  (the  larger  the  N)  covering  the  range  from  Fn).ix  '«>  F  in  are  needed  to  ob¬ 
tain  a  good  answer.  For  example,  we  find  that  in  the  above  case  the  difference 
in  the  energy  deposition  rate  for  N  12  and  N  =  20  is  -0.25  percent  at  100  km, 

+0.  8  percent  at  110  km  and  +1.  (>  percent  at  104  km. 

11.  COMPARISONS 
1 1 . 1  Preliminary  Remark* 

In  this  section,  we  compare  our  analytical  results  to  other  transport  calcula¬ 
tions.  Quantities  of  interest  are  the  electron  flux,  its  hem ispherically  averaged 
value,  and  the  energy  deposition  rate.  Selected  comparisons  will  be  made  for 
monoenergetic  and  Maxwellian  sources. 

The  particular  models  we  use  for  comparisons  are  the  llees.  Range,  and 
Strickland  models.  A  Range  model  was  specifically  developed  for  this  wor  k  and 
to  follow  will  be  a  brief  description  of  its  contents.  We  start  with  the  loss 
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function,  L(E),  which  one  can  obtain  from  the  Bethe  formula  for  energies  above 
several  hundred  eV  or  from  cross  sections  using 


ME)  =  £  WkQexk(E)  +  Z  -  lOoj^iE,  E')  .  C<1) 

k  l 

ME)  has  been  calculated  from  the  cross  sections  in  Strickland  et  al,  for  N.j,  and 
from  current  tabulations  for  0,,  and  0.  *  The  formula  needed  to  specify  the 
energy  loss  is 

7 

o 

Mz.Ko.Mo>  Eo-£  /  d7'4(-11I.1U':M,K0,M0))njM)  ,  (72) 

j  7 

where  is  the  cosine  of  incident  angle,  z  is  aimude,  E  is  the  starting  energy 
at  zt>,  and  n.  (z)  is  the  particle  density  of  the  jth  neutral  species.  The  energy 
deposition  rate  is 


n  {,,(z) 


/  dE  / 

W  O  » 


dn 


I 

,i 


n.(z)  4>(Eo, 


^V*7-  >V 


(7:i) 


For  all  four  models  compared  in  this  section  the  Jacchia  (1D00U)“*  mode) 
atmosphere  was  used.  For  the  Strickland  and  llange  models  the  same  individual 
inelastic  cross  sections  and  constituent  neutral  densities  were  used  and  for  the 
analytic  model  the  same  total  inelastic  cross  section  was  used.  The  weighted 
total  inelastic  cross  section,  Q(K),  was  2.33,  1.31,  0.  atiS,  0.221,  0.0893,  and 
0.0297  times  10  in  cm"  for  0.3,  1,  :i,  10,  30,  and  100  KeV,  respectively. 

The  value  used  for  E^  was  34  eV  and  the  values  used  for  \Y  are  given  in  Section  4. 


11.2  Comparisons  for  Isotropir  iMonoenorgetic  incident  Fluxes 

The  models  to  be  considered  are  the  analytic,  Kees,  and  llange  models. 

Incident  energies  to  be  considered  are  2,  5,  and  10  KeV.  Figure  5  shows  17. ,  vs  7 

•>  '■ 

for  incident  power  densities  of  1  erg/cm“s.  The  Rees  results  were  generated  by 
us  using  his  energy  dissipation  function.  The  analytic  and  llange  results  are  sim¬ 
ilar  as  expected  since  the  corresponding  models  both  contain  the  forward 


23.  Jacchia,  L.G.  (1977)  Thermospheric  temperature,  density,  and  compos. t ion: 
new  models.  Rep.  375,  Smithson.  Astrophys.  Observ. ,  Cambridge,  Mass. 


(kml  l  (km,  Z  (km) 


scattering  approximation.  Departures  result  from  the  differences  in  the  descrip 
tii>n  of  energy  loss.  The  analytic  model  allows  fo  discrete  energy  loss  and  in 
turn  straggling.  Greater  differences  exist  between  the  analytic  and  Rees  results 
the  latter  showing  a  broader  deposition  prof  with  less  total  energy  deposition. 
We  expect  these  latter  r  esults  to  be  the  more  accurate  sinre  they  include  the 
effect  of  pitch  tingle  scattering. 

Il.lt  ('.oniparisoiis  for  Isotropic M;i\ w r Ilian  Incident  h  luvs 

The  models  to  he  considered  are  the  analytic.  Range,  and  Strickland  models 
Characteristic  energies  to  be  considered  tire  1,  2,  and  f>  KcV.  Figure  li  shows 

the  incident  isotropic-Muxwcllian  distributions  for  these  energies.  .Again,  the 

2 

incident  power  densities  are  1  erg  rm's,  Figure  7  shows  rj|,  vs  z  for  the  three 
models  and  three  incident  fluxes.  These  results  exhibit  basically  the  same 
behavior  as  those  appearing  in  Figure  5.  The  Strickland  r  esults  give  a  broader 
distribution  with  a  smaller  total  energy  deposi  ion  rate  due  to  pitch  angle  scatter¬ 
ing,  contained  in  his  model.  For  all  the  analytical  results  of  this  section,  12 
pseudoscatterings  were  used  spanning  the  range  from  1  to  13  KcV.  Based  upon 
our  studies  12  pseudoscatterings  give  a  solution  >f  Eq.  (lf>)  to  within  about 
2  percent  for  all  altitudes  above  104  km. 


ENERGY  (KeV) 


Figure  6.  Incident  Maxwellian  Energy  Distributions 
for  Characteristic  Energies  of  5,  2,  and  1  KeV  Each 
Containing  1  erg/cm2s 
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Figure  7.  Energy  Deposition  Rates  from  the 
Analytic,  Range,  and  Strickland  Models  for 
the  5,  2,  and  1  KeV  Isotropic -Maxwellian 
Sources  as  in  Figure  (> 
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We  are  interested  in  knowing  how  well  the  analytic  model  does  on  the  electron 
flux  itself.  Since  the  model  contains  the  forward  scattering  approximation,  we 
are  specifically  referring  to  ‘he  downward  moving  flux.  To  investigate  this  point, 
fluxes  from  the  analytic  and  Strickland  models  have  been  compared  al  selected 
altitudes  and  energies.  Figure  8  allows  for  such  a  comparison  at  110  km,  for 
the  2  KeV  isotropic-Maxwellian  case.  As  expected,  the  analytic  results  are  high 
near  p  ~  1  md  low  as  «  ->  0  due  to  the  forward-scattering  approximation.  The 
differences  increase  with  decreasing  energy  since  scattering  becomes  more  im¬ 
portant  at  low  energy.  Figure  9  shows  the  hemispherically  averaged  flux  for  the 
two  models  from  ~  1  to  10  KeV  at  110  and  IfiO  km. 


12.  CONCLUSION 

From  the  preceding  results  we  conclude  that  the  forward  scattering  and  aver¬ 
age  discrete  energy-loss  approximations,  which  at  first  seem  quite  severe,  to¬ 
gether  with  the  pseudopartiele  method  of  approximating  the  sums,  produce  a  sol¬ 
ution  for  the  primary  auroral  electron  flux  in  the  downward  hemisphere  which  is 
quite  simple  to  use  and  is  surprisingly  close  to  most  of  the  detailed  numerical 
results  of  Strickland  et  al.  12  The  formulae  we  give  for  the  primary  auroral 
electron  flux  and  the  quantities  derived  from  it  give  us  insight  into  the  physics  of 
the  electron  precipitation  process  and  also  provide  a  means  of  making  rapid  cal¬ 
culations  for  the  purpose  of  analyzing  auroral  data.  We  also  hope  that  the  reader 
will  find  he  results  presented  in  this  paper  useful  in  further  auroral  studies. 
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E  • 10  K«V 


Figure  8.  Differential  Electron  Fluxes  from  the  Analytic  (N  ■  12)  and  Strickland 
Models  at  110  km  for  the  2  KeV  Isotropic -Maxwellian  Source 
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2  KeV  MAXWELLIAN 


Figure  9.  Ileniisphorically  Averaged  Fluxes 
from  the  Analytic  (N  =  12)  and  Strickland 
Models  at  110  and  lf>0  km  for  the  2  KeV  Iso¬ 
tropic -Maxwellian  Source 
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